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By the use of extensive numerical simulations we show that the nearest-neighbor energy level 
spacing distribution P(s) and the entropic eigenfunction localization length of the adjacency matrices 
of Erdos-Renyi (ER) fully random networks are universal for fixed average degree £ = aN ( a and 
N being the average network connectivity and the network size, respectively). We also demonstrate 
that Brody distribution characterizes well P(s) in the transition from a = 0, when the vertices in 
the network are isolated, to a = 1, when the network is fully connected. Moreover, we explore 
the validity of our findings when relaxing the randomness of our network model and show that, 
in contrast to standard ER networks, ER networks with diagonal disorder also show universality. 

Finally, we also discuss the spectral and eigenfunction properties of small-world networks. 

PACS numbers: 64.60.Aq 89.75.Da 


I. INTRODUCTION 


Networks have been used to represent the organiza¬ 
tion of complex systems, such as social networks, Inter¬ 
net, and ecosystems 00 - Depending on the application, 
vertices and edges have different meanings 0,0. For ex¬ 
ample, in condensed matter physics, vertices and edges of 
the ordered network known as Anderson’s tight-binding 
model are the sites and hopping integrals, respectively 
4], Networks can be deterministic, fractal, or random 
5(. Deterministic and fractal networks are constructed 
following specific rules, whereas for random networks a 
set of parameters take fixed values but the network itself 
has a random organization. In this later case it is mean¬ 
ingless to study a single random network, instead a sta¬ 
tistical analysis of an ensemble of networks with the same 
average properties should be performed. Several models 
of random networks have been introduced 0 0 , includ¬ 
ing Erdos-Renyi (ER) random graphs, the scale-free net¬ 
work model of Barabasi and Albert, and the small-world 
networks of Watts and Strogatz. These models are con¬ 
sidered to reproduce the organization of real-world net¬ 
works, such as the Internet, power-grids, and social and 
biological networks 010. Although random graphs 
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fail in predicting most of the properties observed in real- 
world networks, such as power-law degree distributions 
and nonvanishing clustering coefficient 0, such graphs 
have been deeply studied theoretically (e.g. Q). Indeed, 
many results, such as the emergence of percolation, can 
be obtained analytically in ER networks 0, 0 • 


Here we consider the ER random graph model, which 
was introduced by Solomonoff and Rapoport 0 and 
deeply studied later by Erdos and Renyi |, E3. This 
model is also known as uncorrelated random graph 
model. ER networks are constructed by starting with 
N isolated vertices and afterwards, each pair of vertices 
is connected according to a probability a. This process 
is a type of A r2 -realizations Bernoulli process with proba¬ 
bility of success a. Therefore, the number of connections 
follows a binomial distribution. Nevertheless, most real¬ 
izations of this model take into account large values of 
N and small values of a. In this way, the degree distri¬ 
bution tends to a Poisson distribution due to the law of 
rare events. 


Independently of the field, classification, or applica¬ 
tion, a commonly accepted mathematical representation 
of a network is the adjacency matrix. The adjacency 
matrix A of a simple network, i.e. a network having no 
multiple edges or self-edges, is the matrix with elements 
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Ajj defined as 0| 

. _ J 1 if there is an edge between vertices i and j. 

~ 0 otherwise. 

( 1 ) 

This prescription produces TV x N symmetric sparse ma¬ 
trices with zero diagonal elements, where N is the num¬ 
ber of vertices of the corresponding network. The spar¬ 
sity of A is quantified by the parameter a which is the 
fraction of non-vanishing off-diagonal adjacency matrix 
elements. Vertices are isolated when a = 0, whereas 
the network is fully connected for a = 1. Once the ad¬ 
jacency matrix of a network is constructed, it is quite 
natural to ask about its spectral and eigenfunction prop¬ 
erties, which is the main subject of this paper. As com¬ 
monly used, we refer to the spectral and eigenfunction 
properties of the adjacency matrix as the spectral and 
eigenfunction properties of the respective network. 

Moreover, there is a one-to-one correspondence be¬ 
tween the adjacency matrix A and the Hamiltonian ma¬ 
trix H of a ^-dimensional solid, described by Ander¬ 
son’s tight-binding model jj] with zero on-site potentials 
(Hu = 0) and constant hopping integrals ( H\j = 1). Here 
£ is proportional to the average non-zero off-diagonal ad¬ 
jacency matrix elements per matrix row and, therefore, 
may be regarded as the effective dimension of the network 
represented by A; as discussed in [ll| from a random ma¬ 
trix theory (RMT) point of view. This correspondence 
enables the direct application of studies originally moti¬ 
vated on physical systems, represented by Hamiltonian 
sparse random matrices, to complex networks. In the 
network literature, £ is known as the average degree of a 
network. Moreover, notice that 

i = « x N , (2) 

where £ was defined as the mean number of nonzero el¬ 
ements per matrix row. From a mathematical-physicist 
point of view, in the frame of RMT, in 1988 Rodgers and 
Bray |T~2| proposed an ensemble of sparse random ma¬ 
trices characterized by the connectivity £. Since then, 
several papers have been devoted to analytical and nu¬ 
merical studies of sparse symmetric random matrices (see 
for example [Tll4251 j. 

Among the most relevant results of these studies we 
can mention that: (i) in the very sparse limit, £ —> 1, 
the density of states was found to deviate from the 
Wigner semicircle law with the appearance of singu¬ 
larities, around and at the band center, and tails be¬ 
yond the semicircle fl2| - l2Tl |: (ii) a delocalization tran¬ 
sition was found at £ ~ 1.4 [fH - IrI |22| |; (iii) the 
nearest-neighbor energy level spacing distribution P(s) 
was found to evolve from the Poisson to the Gaussian 
Orthogonal Ensemble (GOE) predictions for increasing 
£ [H!, Q, (3 (the same transition was reported for the 
number variance in Ref. JTT( ). More recently, the first 
eigenvalue/eigenfunction problem was also addressed in 
Ref. [23}. Also, non-Hermitian sparse matrices were ap¬ 
proached in Refs. j24j. 


It is relevant to emphasize that the RMT model of 
sparse matrices introduced by Rodgers and Bray [l2j is 
equivalent to adjacency matrices of ER-type networks. 
In fact, motivated by this equivalency and based on 
the ER model, here we study spectral and eigenfunction 
properties of the following random network model: Start¬ 
ing with the standard ER network, we add to it self-edges 
and further consider all edges to have random strengths 
[26j . We call this model as the ER fully random net¬ 
work model. The sparsity a is defined as the fraction of 
the N(N — l)/2 independent non-vanishing off-diagonal 
adjacency matrix elements. Then, as in the standard 
ER model, the ER fully random network model is char¬ 
acterized only by the parameters N and a. However, 
the corresponding adjacency matrices come from the en¬ 
semble of N x N sparse real symmetric matrices whose 
non-vanishing elements are statistically independent ran¬ 
dom variables drawn from a normal distribution with 
zero mean {Aij) = 0 and variance (|A^j 2 ) = (1 + <%)/2. 
According to this definition a diagonal random matrix is 
obtained for a = 0 (Poisson case), whereas the GOE is 
recovered when a = 1. 

Our motivation to study spectral and eigenfunction 
properties of the ER fully random network model is 
twofold. On the one hand, Jackson et al. EH showed 
that the disorder parameter £ fixes some spectral fea¬ 
tures of sparse random matrices equivalent to the adja¬ 
cency matrices we study here. Moreover, it is also known 
that the average degree £ fixes some properties of random 
graphs [27]. On the other hand, there is a large number 
of papers studying spectral and eigenfunction properties 
of complex networks [2Sl [tt! . However, most of those 
studies focus on networks with specific combinations of 
N and a. For instance Refs. [331 - 135 . Hoj investigated 
random networks using concepts of RMT centering the 
attention on networks with fixed N and a, even when 
more sophisticated topological properties are included in 
the network models, such as non-vanishing clustering co¬ 
efficient HH and modular structure [37J. Moreover, Palla 
and Vallay [H| investigated spectral properties of stan¬ 
dard ER networks close to the critical point of percolation 
restricting the analysis to values of 1 < £ < 2.4. 

Then, we realized that the universality (in the sense of 
identifying relevant network parameters) of the spectral 
and eigenfunction properties of ER-type random net¬ 
works has not been completely explored yet. Thus, in 
this paper we undertake this task having as a reference 
RMT models and predictions. Furthermore, our inves¬ 
tigations generalize previous results in the literature by 
not restricting our analysis to a specific combination of 
network parameters or even to a specific network model. 

In the next section we analyze P(s), i.e. the nearest- 
neighbor energy level spacing distribution, the average 
Shannon entropy (S), and the entropic eigenfunction lo¬ 
calization length for the ER fully random network 
model as a function of the connectivity a. It is well 
known that P(s), the probability distribution function of 
the spacings between adjacent eigenvalues, plays a promi- 
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nent role in the description of disordered and quantized 
classically-chaotic systems [HU, [59j. We show that P(s), 
(S), and £n are all invariant for fixed average degree £. 
In addition, we show that Brody distribution fits well 
P(s) in the transition from isolated to fully connected 
networks. By noticing that the ER fully random network 
model displays maximal disorder, in Sect. Ill we explore 
the validity of our findings when relaxing the random¬ 
ness of this network model. Therefore we show that the 
ER network model with diagonal disorder (i.e. ER model 
including random strength self-edges) also exhibits uni¬ 
versality for fixed £. We also apply our approach to a dif¬ 
ferent random network model, so in Sect. IV we comment 
on the spectral and eigenfunction properties of small- 
world random networks. Our conclusions are developed 
in Sect. V. 

II. ERDOS-RENYI FULLY RANDOM 
NETWORKS 

In the following we use exact numerical diagonaliza- 
tion to obtain the eigenvalues E m and eigenfunctions 4' m 
(m = 1 ... N) of the adjacency matrices of large ensem¬ 
bles of random networks characterized by N and a. 

A. Nearest-neighbor energy level spacing 
distribution 

Figure [T] presents the nearest-neighbor energy level 
spacing distribution P(s) for the adjacency matrices of 
ER fully random networks of size N = 1000 and dif¬ 
ferent connectivity values a. The histograms were con¬ 
structed by the use of the 500 unfolded spacings i58j . 
s m = ( E m+1 — E m )/ A, around the band center of 10^ 
random matrices. Here, A is the mean level spacing com¬ 
puted for each adjacency matrix as the slope of the curve 
E m vs. m around the band center. When other network 
sizes are considered, we always construct P(s) from half 
of the total eigenvalues around the band center, where 
the density of states is approximately constant. 

For a = 0, i.e. when the vertices in the network are iso¬ 
lated, the corresponding adjacency matrices are diagonal 
and P(s) follows the exponential distribution, 

P(s) = exp(-s) , (3) 

better known in RMT as Poisson distribution or the spac¬ 
ing rule for random levels [HU. In the opposite limit, 
a = 1, when the network is fully connected, the adja¬ 
cency matrices become members of the GOE (full real 
symmetric random matrices ) an d P(s) closely follows the 
Wigner-Dyson distribution [HI, [&Qj, 

p {s) = ^sexp(-^s 2 ) . (4) 

Then, by increasing a from zero to one, the shape of 
P(s) should evolve from the Poisson to the Wigner-Dyson 


distributions. This transition, partially observed for our 
random network model in Refs. EMI from a pure RMT 
point of view, is well depicted in Fig.Q] where we also plot 
Eqs. ([3]) and (Q} as reference. In fact, it is interesting to 
mention that for relatively small values of a the limiting 
GOE statistics is recovered. The transition from Poisson 
to Wigner-Dyson in the spectral statistics has also been 
reported for adjacency matrices corresponding to other 
complex network models in Refs. [30M37I1 . 

Here, in order to characterize the shape of P(s) for 
our random networks we use the Brody distribution 
f6H . [62 ] which was originally derived to provide an in¬ 
terpolation expression for P(s) in the transition from 
Poisson to Wigner-Dyson distributions by making the 
ansatz P(s) = Cis^ exp (—C 2 S /3+1 ) (with Ci^ depending 
on ft). Then, after proper normalization Brody distribu¬ 
tion reads [Hi], IH^] 

P(a) = {P + 1 )a/?s /3 exp (~aps p+1 ) , (5) 

where ag = [P(/3 + 2/ft + 1)]^ +1 , T(-) is the gamma func¬ 
tion, and ft, known as Brody parameter, takes values in 
the range [0,1]. ft = 0 and /? = 1 reproduce the Pois¬ 
son and Wigner-Dyson distributions, respectively. We 
want to remark that, even though Brody distribution 
has been extensively used to characterize P(s) having 
fractional power-law level repulsion [HU, it has been ob¬ 
tained through a purely phenomenological approach and 
the Brody parameter has no decisive physical meaning 
[64j but serves as a measure for the degree of mixing be¬ 
tween the Poisson and GOE statistics. In particular, as 
we show below, the Brody parameter will allows us to 
identify the onset of the delocalization transition and the 
onset of the GOE limit in our random network models. 

Black dashed lines in Fig. |T] represent the fittings of 
the numerically obtained P(s) with Eq. (fo[). The fitted 
values of ft are given in the corresponding panels. This 
figure shows that Brody distribution provides very good 
fittings for the P{s) of the adjacency matrices of ER fully 
random networks. In fact, Brody distribution also works 
well for other complex networks models IHHl fibl . 

Now we also construct histograms of P(s) for a large 
number of values of a to extract systematically the corre¬ 
sponding values of ft. Figure Ufa) reports /3 versus a for 
five different network sizes. Notice that in all five cases 
the behavior of ft is similar: ft shows a smooth transi¬ 
tion from zero (Poisson regime) to one (Wigner-Dyson 
or GOE regime) when a increases from a -C 1 (mostly 
isolated vertices) to one (fully connected networks). No¬ 
tice also that the larger the network size IV, the smaller 
the value of a needed to approach the GOE limit. This 
is, in fact, the reason we observed in Fig. |T] that P(s) is 
very close to the Wigner-Dyson distribution already for 
a = 0.00426 < 1. 

We recall that the parameter £ (see Eq. Q) fixes the 
density of states of sparse random matrices [ll] equiva¬ 
lent to the density of connections in the adjacency matri¬ 
ces we study here. Moreover, £ was defined as the effec¬ 
tive dimension of the sparse random matrix model E2- 
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FIG. 1: (Color online) Nearest-neighbor energy level spacing 
distribution P(s) for ER fully random networks of size N = 
1000 and different connectivity values a (red histograms). 
Blue lines correspond to Poisson and Wigner-Dyson distribu¬ 
tion functions given by Eqs. @ and 0, respectively. Black 
dashed lines are fittings of the histograms with the Brody 
distribution of Eq. ©, where the fitted values of fj are given 
in the corresponding panels. The histograms were computed 
from 5 x 10 5 unfolded spacings. 

Also, Martfncz-Mendoza et al. (65| showed that scatter¬ 
ing and transport properties of tight-binding ER fully 
random networks are universal for fixed disorder param¬ 
eter £. Thus, it make sense to explore the dependence of 
f3 on the average degree £; in fact, in Fig. [2])b) we show 
this dependency. We observe that curves for different 
network sizes N fall on top of a universal curve. This 
means that once £ is fixed, no matter the network size, 
the shape of P(s ) is also fixed. We also note that the 
transition in the form of P(s) takes place in the interval 
1 < £ < 7; i.e. when £ < 1 (£ > 7), P(s) has the Poisson 
(Wigner-Dyson) shape. 

Therefore, we verify the invariance of the form of P(s) 
for fixed average degree £ by (i) choosing four represen¬ 
tative values of /? (0.25, 0.5, 0.75, and 0.98 ~ 1); (ii) 
extracting the corresponding values of £ from the univer¬ 
sal curve of Fig. ©b); and (iii) constructing histograms 
of P(s) for several network sizes for each of the chosen 
values of £. Figure [3] shows that once £ is fixed, the form 
of P(s) is invariant. In addition, we include in each panel 
the corresponding Brody distributions. 

Notice that Fig.©b) also provides a way to predict the 
shape of P(s) of ER fully random networks once the aver¬ 
age degree £ is known: When £ < 1, P(s) has the Poisson 
shape. For £ > 7, P(s) is practically given by the Wigner- 
Dyson distribution. While in the regime 1 < £ < 7, P(s) 
is well described by Brody distributions characterized by 
a value of 0 < /3 < 1. Thus, £ = 1 and £ = 7 mark the 
onset of the delocalization transition and the onset of the 
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FIG. 2: Brody parameter f) as a function of (a) the connec¬ 
tivity a and (b) the average degree £ = aN for ER fully 
random networks of sizes ranging from N = 250 to 4000. 
Dashed vertical lines at £ = 1 and £ = 7 mark the onset of 
the delocalization transition and the onset of the GOE limit, 
respectively. Error bars are not shown since they are much 
smaller than symbol size. 

GOE limit, respectively. 

B. Entropic eigenfunction localization length 

In order to characterize quantitatively the complex¬ 
ity of the eigenfunctions of random matrices (and of 
Hamiltonians corresponding to disordered and quantized 
chaotic systems) two quantities are mostly used: (i) the 
information or Shannon entropy and (ii) the eigenfunc¬ 
tion participation number. These measures provide the 
number of principal components of an eigenfunction in a 
given basis. In fact, both quantities have been already 
used to characterize the eigenfunctions of the adjacency 
matrices of random network models (see some examples 

in Refs. m he, in in a an aa ). 

Here, we use the Shannon entropy, which for the eigen¬ 
function 'f> m is given as 

N 

S=-J2(*n) 2 H*n) 2 ■ (6) 

n —1 

S allows us to compute the so called entropic eigenfunc¬ 
tion localization length [7Gj, i.e. 

l N = A^exp [- (S GOE - (S))] , (7) 

where S GOE ~ ln(IV/2.07) is the entropy of a random 
eigenfunction with Gaussian distributed amplitudes. We 
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FIG. 3: (Color online) Nearest-neighbor energy level spac¬ 
ing distribution P(s) for ER fully random networks of sizes 
ranging from N = 250 to 4000 and different values of the av¬ 
erage degree £ (color histograms). Black dashed lines are the 
Brody distribution of Eq. (0 with /3 = 0.25, 0.5, 0.75, and 
0.98. Each histogram was computed from 5 x 10 5 unfolded 
spacings. 



FIG. 4: Average Shannon entropy (S) normalized to Sgoe 
as a function of (a) the connectivity a and (b) the average 
degree £ for ER fully random networks of sizes ranging from 
N = 500 to 4000. Each point was computed by averaging 
over 10 6 eigenfunctions. 


average over all eigenfunctions of an ensemble of adja¬ 
cency matrices of size N to compute (S) With this 
definition, when a = 0, since the eigenfunctions of the 
adjacency matrices of our random network model have 
only one non-vanishing component with magnitude equal 
to one, ( S) = 0 and In ~ 2.07. On the other hand, for 
a = 1, ( S) = Sgoe and the fully chaotic eigenfunctions 
extend over the N available vertices in the network, i.e. 
£ n ~N. 

Figures 0a) and 0(a) show (S) / S GOE and £n/N, re¬ 
spectively, as a function of the connectivity a for the 
adjacency matrices of ER fully random networks of sizes 
N = 500, 1000, 2000, and 4000. We observe that the 
curves (S) /S GO e and £n/N have the same functional 
form as a function of a. Notice that this behavior is also 
observed for /? in Fig. 0a). Also these curves are dis¬ 
placed to the left for increasing N. On the other hand, 
when we plot (S) /S GOE and In /N as a function of the 
average degree £, we observe the coalescence of all curves 
into universal ones (see Figs. 0b) and 0b)). In addi¬ 
tion, note that the point at which every ER fully random 
network becomes globally connected (in the sense that it 
does not contain isolated sub-networks), a = (In N)/N 
[22}, occurs when £n/N « 1/2 (see the dashed vertical 
lines in Fig. 0a)). 

From Fig. 0b) it is clear that the universal behav¬ 
ior of the curve £n/N as a function of the average de¬ 
gree £ can be easily described: (i) £n/N transits from 
s=s 2.07/N ~ 0 to one by moving £ from zero to TV; 
(ii) for £ ~ 2 the eigenfunctions are practically local¬ 


ized since £n ~ 1; hence the delocalization transition 
takes place around £ « 2, which is close to previous the¬ 
oretical and numerical estimations l-B-lla. 1221 : and (iii) 
for £ > 200 the eigenfunctions are practically chaotic 
and fully extended since £n « N. Despite the fact that 
the transition region for £n/N , which takes place in the 
range 2 < £ < 200, is very large as compared to that for 
/3, both quantities are highly correlated and characterize 
well the Poisson to Wigner-Dyson transition of spectral 
and eigenfunction properties of the adjacency matrices of 
ER fully random networks as a function of £. 

III. OTHER ERDOS-RENYI RANDOM 
NETWORKS 

Notice that with the prescription given above our ER 
fully random network model displays maximal disorder, 
because averaging over the network ensemble implies av¬ 
erage over connectivity and over connection strengths. 
With this averaging procedure we get rid off any individ¬ 
ual network characteristic (such as scars l72l| . which in 
turn produce topological resonances 0 ) that may lead 
to deviations from RMT predictions used here as a ref¬ 
erence (see also 0 )- More specifically, we choose this 
network model to retrieve well known random matrices 
in the appropriate limits — remember that a diagonal 
random matrix is obtained for a = 0, when the vertices 
are isolated, whereas a member of the GOE is recovered 
for a = 1, when the network is fully connected. 
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FIG. 5: Entropic eigenfunction localization length In nor¬ 
malized to IV as a function of (a) the connectivity a and (b) 
the average degree £ for ER fully random networks of sizes 
ranging from N = 500 to 4000. Dashed vertical lines in (a) 
mark the values of a = (In N)/N at which e very ER fully 
random network becomes globally connected [27j]. Dashed 
vertical lines in (b) at £ = 2 and £ = 200 mark the onset of 
the delocalization transition and the onset of the GOE limit, 
respectively. Each point was computed by averaging over 10 6 
eigenfunctions. 


However, it is important to add that the maxi¬ 
mal disorder we consider above is not necessary for a 
graph/network to exhibit universal RMT behavior. In 
fact, it is well known that tight-binding cubic lattices 
with on-site disorder (known as the three-dimensional 
Anderson model i)> forming networks with fixed reg¬ 
ular connectivity having very dilute adjacency matrices, 
show RMT behavior in the metallic phase (see for exam¬ 
ple Refs. [zilzi). Moreover, it has been demonstrated 
numerically and theoretically that graphs with fixed con¬ 
nectivity show spectral [z3 III and scattering [FI, [§d] 
universal properties corresponding to RMT predictions. 
In this case the disorder is introduced either by choosing 
random bond lengths [z3, [zl, H (which is a parame¬ 
ter not present in our network model) or by randomizing 
the vertex-scattering matrices [78 j (somehow equivalent 
to consider random connection strengths). Some of the 
RMT properties of quantum graphs have already been 
tested experimentally by the use of small ensembles of 
small microwave networks with fixed connectivity 8lj . 
Furthermore, complex networks having specific topolog¬ 
ical properties (such as small-world and scale-free net¬ 
works, where randomness is applied only to the connec¬ 
tivity) show signatures of RMT behavior in their spectral 
and eigenfunction properties (34M361 HH, HH . 


Therefore, in the following we search for the scaling 
properties, if any, of ER-type random networks when 
the condition of maximal disorder considered above is 
relaxed. 


A. Standard Erdos-Renyi random networks 

In the standard ER random network model I1-E3, 
the corresponding adjacency matrices are random ma¬ 
trices with zeros in the main diagonal and ones as non¬ 
vanishing off-diagonal elements; i.e. in the adjacency ma¬ 
trix vertices and edges are represented with zeros and 
ones, respectively (see Eq. (JT]) ) . Even though it has been 
shown that the P(s) of standard ER random networks is 
close to the Wigner-Dyson shape for large connectivity 
(a —>• 1) [3J, [35|, [38}, notice that P(s ) can not show the 
Poisson to Wigner-Dyson transition since in the limit of 
vanishing connectivity (a = 0) the corresponding adja¬ 
cency matrices are the null matrix. 

Figure [5] shows the distribution P(s) for standard ER 
random networks of size N = 1000. We observe that once 
£ is large enough P(s) acquires the expected Wigner- 
Dyson shape. However, for smaller values of £, P(s) de¬ 
velops two components: (i) a prominent peak at s = 0 
and (ii) a broad part having a local maximum at s > 0. 
The presence of these two components avoids the use of 
the Brody distribution to describe the shape of P(s) in 
the transition from isolated vertices to fully connected 
networks. The same behavior is observed for any other 
value of N. 

It is fair to say that the Brody parameter has already 
been reported as a function of £ for standard ER random 
networks in Ref. [28J; however, for 1 < £ < 2.4 only. 
The reason to report such a narrow range of £ is because 
for vanishing connectivity this model does not approach 
the Poisson limit and P(s) can not be fitted by Brody 
distribution. 

Concerning the Shannon entropy and the entropic 
eigenfunction localization length for this random net¬ 
work model, we observe, as expected, that ( S) /S GOE and 
In /N approach zero when a —> 0, whereas they approach 
one when a —>• 1. Thus, the delocalization transition of 
this network model can be well characterized by both 
(S) and In- However, there is no universal scaling of 
( S) /S'qoe and In/N when plotted as a function of £ (not 
shown here). 

B. Erdos-Renyi random networks with diagonal 
disorder 

We construct the ER random network model with diag¬ 
onal disorder by adding self-edges with random strengths 
such that the corresponding adjacency matrices acquire 
statistically independent random variables (drawn from 
a normal distribution with zero mean and variance one) 
in their main diagonal [82|. With this construction, we 
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FIG. 6: (Color online) Nearest-neighbor energy level spacing 
distribution P(s) for standard ER random networks of size 
N = 1000 and different values of the average degree £ (color 
histograms). Black dashed line is the Wigner-Dyson distribu¬ 
tion of Eq. ©■ The arrow indicates the direction of increasing 
£. Each histogram was computed from 5 x 10 5 unfolded spac- 
ings. 
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do observe a clear Poisson to Wigner-Dyson transition in 
the form of P(s ) when a moves from zero to one. More¬ 
over, P(s) can be effectively fitted by the Brody distri¬ 
bution (not shown here). Thus, Fig. [3(a) presents the 
Brody parameter /3 as a function of the average degree £ 
for ER random networks with diagonal disorder. As well 
as for ER fully random networks, here we observe that 
the curves of /3 vs. £ for different network sizes N fall 
on top of a universal curve. In addition, we also observe 
universal behavior for the entropic eigenfunction local¬ 
ization length normalized to N, as a function of £ 
(see Fig. [3(b)). The delocalization transition in ER net¬ 
works with diagonal disorder has also been investigated 
in Ref. [32|, but as a function of the disorder strength, 
by the use of P(s) for networks having 3 < £ < 10. 

It is relevant to add that even though we observe the 
collapse of the curves (3 vs. £ and £n/N vs. £ for ER ran¬ 
dom networks with diagonal disorder of different sizes, 
these universal curves are slightly displaced to the left 
as compared to the corresponding curves for ER fully 
random networks. More specifically, the onset of the de- 
localization transition and the onset of the GOE limit 
occurs for smaller values of £ in the case of ER random 
networks with diagonal disorder. See dashed lines and 
red curves in Fig. [3 (included for comparison purposes) 
which correspond to ER fully random networks. 

IV. SMALL-WORLD NETWORKS 

Once we analyzed the universal properties of ER-type 
random networks, it make sense to further explore other 
random networks models to look for universal properties. 
For this task we consider small-world (SW) networks InI . 
A SW network, as defined in 1998 by Watts and Strogatz 
[IU, is constructed by randomly rewiring the edges of a 
regular ring network consisting of N vertices connected 
to their kj2 nearest-neighbors (k > 2 must be an even 
number) |8J] . Then, for every vertex, every right-handed 
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FIG. 7: (Color online) (a) Brody parameter j3 and (b) entropic 
eigenfunction localization length In (normalized to N) as a 
function of the average degree £ for ER random networks 
with diagonal disorder of sizes ranging from N = 500 to 4000. 
Dashed red vertical lines, displayed for comparison purposes, 
at £ = 1 and £ = 7 in (a) [at £ = 2 and £ = 200 in (b)] 
mark the onset of the delocalization transition and the onset 
of the GOE limit, respectively, of ER fully random networks; 
see Fig. ©b) [see Fig. EKb)]. Red data sets corresponding to 
ER fully random networks of size N = 4000 are also included 
as a reference. 


edge is reconnected with probability p to a vertex chosen 
uniformly at random. In the standard SW model multi¬ 
ply connected vertex pairs and self-connections are not 
allowed. Notice that for p = 0 the SW network becomes 
the original regular ring; whereas for p = 1, a random 
network is obtained where every vertex has a minimum 
degree of k/2. 

Notice that: (i) The parameter k is in fact equivalent 
to the average degree £, given in ©, which fixes the 
spectral and eigenfunction properties of ER fully ran¬ 
dom networks and ER random networks with diagonal 
disorder. However, k cannot take non-integer values nor 
values less than two as £ does. Also, it is not an average 
quantity, (ii) The rewiring probability p is a parameter, 
independent of fc, which drives the SW network model 
from regular to random. Also, we should stress that as 
well as for standard ER random networks, the adjacency 
matrices of standard SW networks have zeros in their 
main diagonal. This prohibits the use of the Brody dis¬ 
tribution to fit the P(s) of standard SW networks when 
p and k are both small. However, for large p and fc, 
P(s) becomes close to the Wigner-Dyson distribution, as 
shown in Refs. [3dl l36l. l38l. l50j| (for an analytical approach 
to the spectra of SW networks see Ref. [851). 

Then, as we did in the previous Section for ER ran- 















dom networks, here we also consider SW networks with 
diagonal disorder. That is, we include self-edges hav¬ 
ing random strengths (drawn from a normal distribution 
with zero mean and variance one) to each vertex in the 
network, such that the corresponding adjacency matrices 
exhibit random variables in their main diagonal, whereas 
the edges joining vertex pairs are still represented by ones 
in the adjacency matrices. In this way we guarantee that 
P(s) will have the Poisson shape when p and k are both 
small. Moreover, we observe (not shown here) that by 
including diagonal disorder , the Brody distribution fits 
reasonably well the P(s) of SW networks for any combi¬ 
nation of p and k. The delocalization transition in SW 
networks with diagonal disorder has also been investi¬ 
gated in Ref. GOT as a function of p, N , and disorder 
strength by the use of the inverse participation ratio of 
eigenfunctions (equivalent to the entropic eigenfunction 
localization length we use here); however, the goal there 
was the implementation of a quantum algorithm able to 
efficiently simulate the network. 

Figure |5] shows the Brody parameter j3 as a function of 
p (for several values of k ) and k (for several values of p) 
for SW networks with diagonal disorder of size N = 1000. 
Except for k = 2, where P(s) is very close to Poissonian 
for any p (i.e. (3 « 0), since the corresponding adja¬ 
cency matrix is almost diagonal, P(s) shows the transi¬ 
tion from Poisson to Wigner-Dyson as a function of p (see 
Fig. (SKa)). From this figure we can see that the larger 
the value of k, the smaller the value of p needed for /3 to 
approach one. Notice, for example, that for k = 20 the 
value of p = 0.001 -C 1 makes /3 to be very close to one. 
Clearly, the transition from Poisson to Wigner-Dyson in 
the form of P(s ) is also observed as a function of k for 
fixed p (see Fig. [Sfb)). 

Also, in Ref. |34] the Brody parameter has been re¬ 
ported as function of p but for standard SW networks 
with fixed size N = 2000 and average degree k = 40. 
The goal there was to show that /3 is correlated with two 
important network parameters: the characteristic path 
length and the clustering coefficient; quantities which are 
out of the scope of our paper. Anyhow, the delocalization 
transition was clearly shown as a function of p (analogous 
to our curves reported in Fig. Uta) for SW networks with 
diagonal disorder). 

We have to stress that, in contrast to ER random net¬ 
works with diagonal disorder, where the average degree £ 
fixes their spectral and eigenfunction properties, the de¬ 
gree k does not fix the properties of SW networks with 
diagonal disorder. This fact is clearly visible in Fig. [9l 
where we show £jv/ N as a function of p for SW networks 
of four different sizes. Each panel corresponds to a fixed 
value of degree k where the curves approach each other 
only when p > 0.1. 

Finally, we have to mention that in analogy with the 
ER fully random networks we studied in Sect. II, here we 
also considered the case of SW fully random networks; 
i.e. all non-zero entries of the corresponding adjacency 
matrices, including the main diagonal, were considered as 



P 



k 


FIG. 8: Brody parameter /? as a function of (a) the rewiring 
probability p and (b) the degree k for SW networks with di¬ 
agonal disorder of size N = 1000. 



FIG. 9: (Color online) Entropic eigenfunction localization 
length h v (normalized to N) as a function of the rewiring 
probability p for SW networks with diagonal disorder of sizes 
ranging from N = 500 to 4000. (a) k = 2, (b) k = 8, (c) 
k — 12, and (d) k = 20. 


random variables (in fact, this case has been approached 
in Ref. [3(| where the transition from Poisson to Wigner- 
Dyson in the shape of P(s) was reported as a function of 
p for a network with the fixed parameters N = 1600 and 
k = 8). However, we have not observed a substantial 
difference to the results reported in Figs. [5] and [5] (not 
shown here). 
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V. CONCLUSIONS 

We have studied numerically some spectral and eigen¬ 
function properties of Erdos-Renyi-type random net¬ 
works focusing our attention on universality. 

In particular, we have shown for Erdos-Renyi (ER) 
fully random networks (where all non-vanishing adja¬ 
cency matrix elements are Gaussian random variables) 
and ER networks with diagonal disorder (where the di¬ 
agonal adjacency matrix elements are Gaussian random 
variables, whereas the rest of non-vanishing matrix ele¬ 
ments are ones) that: (i) The nearest-neighbor energy 
level spacing distribution P(s), the average Shannon en¬ 
tropy (S'), and the entropic eigenfunction localization 
length are universal for fixed average degree £ = aN 
(where a and N are the network connectivity and the 
network size, respectively); (ii) Brody distribution fits 
well P(s) in the transition from Poisson (a = 0; isolated 
vertices) to Wigner-Dyson (a = 1; fully connected net¬ 
work); and (iii) the Brody parameter /? as a function of 
£ displays an invariant curve. 

This analysis provides a way to predict the shape of 
P(s) of ER-type random networks once the parameter £ 
is known. 

Specifically, for ER fully random networks and ER net¬ 
works with diagonal disorder we have found that when 
0 < £ < 1, P(s) is well described by the Poisson shape. 
This range of £ values coincides with the regime where 
a typical ER random graph is composed of isolated trees 
[271 ] . We have heuristically located the delocalization 


transition point around £ « 1, in close agreement with 
the transition value of £ « 1.4 reported in Refs. 0- 
[lB,EH- Also note that £ « 3.5, known as the average de¬ 
gree value at which the diameter of an ER random graph 
equals the diameter of the giant cluster (27)], is located 
about half the way in the Poisson to Wigner-Dyson tran¬ 
sition. Also, we have observed that for P(s) to approach 
the Wigner-Dyson shape £ > 7 is needed. 

However, we have determined that there is no universal 
scaling of ( S) /S GOE and In/N when plotted as a function 
of £ for the standard ER model nor for SW networks; even 
though these two quantities describe well the delocaliza¬ 
tion transition of both models. Moreover, we can affirm 
that diagonal disorder is needed for a random network 
model to show the Poisson to Wigner-Dyson transition 
in the form of P(s). 
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